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SECTION  1 

DESCRIPTION  OF  1977  ROBIN  PROGRAM 


The  1977  ROBIN  program  is  an  update  of  the  April  1972 
HIROBIN  program.  The  1977  program  contains  several  changes  from 
the  1972  version.  The  changes  were  made  in  order  to  eliminate 
the  Mach  1  and  supersonic  bias,  to  extend  and  modify  the  drag 
table  so  that  the  temperatures  of  sphere  and  Datasonde  flights 
were  in  agreement,  to  improve  the  method  of  calculating  RMS 
wind  and  thermodynamic  errors,  to  improve  the  first  guess  at 
temperature,  and  to  standardize  input  constants.  A  discussion 
of  the  changes  made  and  how  they  were  incorporated  into  the  1977 
ROBIN  program  is  contained  in  the  following  paragraphs. 

Removal  of  Supersonic  and  Mach  1  Bias 

The  1977  ROBIN  program  is  designed  to  remove  both  the  super¬ 
sonic  and  Mach  1  bias  in  all  parameters  -  density,  pressure, 
temperature,  and  winds.  The  technique  used  to  remove  bias  is 
to  simulate  the  bias  by  a  theoretical  trajectory  and  then  subtract 
this  bias  from  the  standard  ROBIN  output.  Thus  the  1977  ROBIN 
(with  respect  to  bias  error  removal)  is  essentially  the  1972 
HIROBIN  with  additional  subroutines  to  calculate  and  subtract 
the  bias  error.  Details  are  as  follows  (see  Figure  1) .  The 
editing  program  and  the  descending  altitude  subroutine  which 
determines  where  data  processing  of  the  flight  begins  have  not 
been  changed  in  the  1977  program.  An  initial  position  and 
velocity  of  the  sphere,  however,  are  calculated  from  the  radar  data 
immediately  after  the  determination  of  five  consecutive  descending 
altitudes.  From  this  initial  position  and  velocity,  a  theoretical 
sphere  trajectory  is  generated  by  assuming  the  sphere  fell  in  the 
1976  Standard  Atmosphere  with  no  winds  present.  The  theoretical 
trajectory  is  terminated  at  60  Km  since  both  Mach  1  and  super¬ 
sonic  bias  occur  above  this  altitude.  The  theoretical  profile 
of  ( t , X, Y, Z )  coordinates  is  then  treated  as  if  it  represented 
radar  coordinates.  Logic  of  the  1972  HIROBIN  program  with  the 


Figure  1.  1977  ROBIN  Program 


19-21,  51-35  filters  is  used  to  generate  atmospheric  parameters 
of  winds,  temperature,  density,  and  pressure  from  the  theoretical 
trajectory.  Any  differences  between  the  computed  parameters  and 
the  1976  Standard  Atmosphere  values  are  the  bias  errors  resulting 
from  the  filters.  These  differences  in  density,  temperature,  and 
winds  are  stored  in  a  bias  correction  array.  Bias  errors  in 
pressure  are  not  stored  since  using  the  ideal  gas  law,  pressure 
can  be  determined  directly  from  temperature  and  density  measure¬ 
ments.  After  generation  of  the  bias  correction  table,  processing 
of  the  radar  data  continues  in  the  same  manner  as  was  done  in 
the  1972  HIROBIN  program.  The  1977  output  is  obtained  by  sub¬ 
tracting  the  bias  correction  table  from  the  radar  processed 
data.  The  correction  is  applied  by  interpolation  of  the  bias 
error  table  at  each  altitude  where  processed  tracking  data 
occurs  -  except  for  altitudes  between  Mach  1.3  and  M'ach  .85. 

During  this  region,  where  the  Mach  1  spike  occurs,  the  bias 
correction  table  is  interpolated  based  on  Mach  number. 

A  feature  of  the  bias  removal  procedure  is  its  ability  to 
remove  both  the  supersonic  bias  and  Mach  1  spike.  Examples 
of  typical  bias  errors  appear  in  Figure  2. 

Drag  Table 

The  drag  table  used  in  the  1977  ROBIN  program  is  an 
extended  and  modified  version  of  the  Bailey  and  Hiatt  drag  table 
used  in  the  1972  HIROBIN.  The  drag  table  was  extended  in  the  region 
adjacent  to  Mach  2  and  Reynolds  Number  1600  so  that  those  high 
apogee  flights  that  in  the  past  went  out  of  the  drag  table  would 
remain  within  given  drag  values.  These  extended  drag  coefficients 

were  obtained  from  Bailey  and  Hiatt's  ballistic  sphere  measure- 

.  [1] 
ments 

Changes  were  made  to  the  subsonic  drag  coefficients 
between  Mach  0  and  0.5  and  Reynolds  Numbers  between  5000  and 
35,000.  The  purpose  of  these  changes,  as  explained  in  Reference  2, 
was  to  bring  sphere  temperature  measurements  into  agreement  with 
Datasonde  measurements  in  the  region  between  32  and  58  Km. 
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The  drag  table  used  in  the  1977  ROBIN  program  is  shown  in 
Figure  3. 

1976  Standard  Atmosphere 

The  calculation  of  pressure  and  temperature  in  the  ROBIN 
program  requires  a  first  guess  of  temperature  (usually  near  95  Km) . 
In  the  1972  HIROBIN  program  the  initial  temperature  was  obtained 
from  the  1962  Standard  Atmosphere.  In  the  1977  ROBIN,  the  first 
guess  is  taken  from  the  1976  Standard  Atmosphere.  The  1976 
Standard  Atmosphere  is  also  used  in  the  1977  ROBIN  program  for 
the  generation  of  the  theoretical  trajectory  from  which  the 
supersonic  and  Mach  1  bias  correction  table  is  generated. 

RMS  Errors 

The  equation  used  in  the  1972  HIROBIN  program  to  estimate 
the  radar  error  failed  to  take  into  account  the  correlation  in 
velocities  introduced  by  the  filtering  process.  Consequently, 
the  error  estimates  were  too  small  by  fifty  percent  or  more  at 
the  higher  altitudes.  Reference  3  discusses  the  correlated 
error  problem  and  provides  an  equation  for  estimating  the  influence 
of  the  correlated  errors.  The  1977  ROBIN  has  been  modified  to 
include  the  correlated  error  term.  Thus,  assuming  that  the 
FPS/16  radar  performs  according  to  specifications,  the  1977 
ROBIN  will  provide  an  accurate  estimate  of  wind  and  thermodynamic 
errors  arising  from  the  radar  track. 

Gravitational  Constant 

To  maintain  consistency  with  other  programs  used  in  the 
reduction  of  meteorological  data,  the  equation  used  to  calculate 
the  gravitational  constant  in  the  1977  ROBIN  program  is 

Gs  =  -9.780356  ^1  +  5.2885  X  10-5  (sin2  (<J>)  )  -  5.9  X  10~8  (sin2  ( 2<J> ) ))  , 

where  v  is  latitude.  No  discernable  change  in  thermodynamic 
or  wind  data  output  results  from  the  modification  of  the  gravi¬ 
tational  constant. 


5 


i H6i  onvi>  3luvi  yvao 


• 

ir 

c 

O  rv 

X 

m 

o- 

ro 

n. 

X 

X 

O 

ro  no  o 

ro 

o  X 

o 

ro 

O' 

n- 

4- 

ro 

c 

— 

o 

CD 

ro 

ro 

ro 

44 

o 

44 

ro 

o 

© 

4- 

ro 

rv 

44 

© 

-r 

o 

r- 

X 

X 

as 

O' 

o 

— i 

ro 

no 

—4  44  —4 

o 

o 

o 

O' 

CD 

r~ 

r- 

ro 

ro 

rv 

—* 

O' 

ro 

— 1 

f- 

ro 

ro 

44 

O' 

O' 

© 

r- 

ro 

ro 

ro 

c 

c 

c 

o 

c 

C 

o 

— 

— 

— 

— 

— 

—  — 

— 

— 

o 

o 

c 

o 

C 

c 

c 

c  o 

O' 

ro 

ro 

X 

X 

X 

X 

X 

ro 

ro 

ro 

«y 

< 

ro 

- 

~ 

- 

- 

~ 

- 

- 

- 

- 

- 

-* 

- 

- 

- 

- 

- 

- 

- 

- 

4 

- 

- 

- 

- 

- 

- 

o 

O 

c 

o  o 

o 

©  o 

o 

o 

O 

o 

o 

O 

o 

o 

O 

o 

C 

o 

• 

X 

o 

ro  O 

ro 

n- 

-T  rJ  —« 

n- 

ro 

o  n> 

CD 

a? 

ao 

o 

ro 

ro 

rv 

no 

X 

CD 

44 

r- 

X 

rv 

X 

X 

o  n- 

ro 

o 

ro 

ro 

ro 

o 

o 

c 

o 

o 

r— 

ro 

ro 

X 

XXX 

■j- 

r* 

ro 

no 

4—4 

C  O'  CD 

ro 

X 

rv 

—4 

O' 

r- 

ro 

O 

rv  r- 

no 

o  n- 

X 

ro 

no 

— 

4- 

©  O- 

© 

f- 

ro 

o 

-* 

*r 

— * 

— 

-* 

4-.  —  -4 

-• 

-* 

— 

-* 

o  o  o 

o 

O  O' 

O' 

o- 

n- 

r*- 

ro 

ro 

ro  x 

X 

X 

X 

X 

X 

ir. 

ro 

ro 

ro 

•» 

4 

- 

- 

- 

- 

- 

~ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

J.-- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

— 

— 

- 

o  o  o  o 

o  o  o 

o 

©  o 

O  O 

c 

o 

C 

o  © 

o 

© 

rv 

* 

e 

a • 

O  ro 

ro 

no 

O  O' 

ro  — 

O' 

ro  r~ 

O 

ro 

ro 

ro  rv  n~ 

X 

C  no 

no  ro 

n- 

n- 

c 

o 

©  r- 

© 

nj 

rv 

©  © 

cc 

X 

n. 

c 

ro 

ro 

ro 

ro 

vT 

■X 

r- 

o  o 

O'  o-  o 

r~ 

n- 

ro 

ro 

X 

ro 

ro 

ro  O  as 

f- 

X 

m 

CD 

rv 

r— 

—4 

r- 

ro 

44 

© 

f'- 

ro 

X 

ro 

ro 

nj 

-4 

o  © 

c 

— 

— 

— 

— 

— 

-« 

— 

ro 

ro 

ro 

ro 

ro 

— 

-  — 

•- 

— 

— ■ 

— 

-* 

o 

C 

O  O  o  O' 

as 

r- 

ro 

ro 

ro 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ro 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

-  -- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

O 

o 

c 

©  o  © 

o 

o 

C 

o 

o 

c 

O 

O  © 

o 

O 

no 

r- 

C  Is- 

ro 

ro 

r~ 

ro 

no 

O  X  — • 

O  x  r- 

o  ro 

ro 

ro 

O  no 

X 

X 

ro  r- 

CD 

r~ 

n- 

n- 

ro 

no 

ro 

® 

© 

O 

ro 

,4 

O  ro 

X 

o 

c 

ro 

ro 

s? 

in 

m 

X)  O' 

X 

ro 

no 

no 

no 

ro 

4-4  4-4  4-< 

4- 

o  o 

O' 

cd  r- 

ro  x 

44 

O' 

r- 

X 

rv 

O  X 

co 

rv  O' 

X. 

rv 

44 

O' 

© 

© 

r- 

ro 

X 

ro 

rv 

ro 

o 

— * 

-* 

-* 

-* 

~4 

-* 

*-• 

rvj 

ro 

no 

ro 

ro 

no 

no 

no 

no 

ro 

ro  ro  nJ 

ro 

ro 

-• 

44 

-* 

o 

o  o 

O  CO 

r~ 

in- 

ro 

ro 

o 

ro 

X 

X 

X 

X 

X 

X 

X 

X 

X 

-* 

- 

- 

- 

- 

- 

~ 

- 

- 

- 

- 

— 

- 

- 

- 

-4 

-4 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

— 

— 

- 

— 

- 

o 

o 

© 

o 

o 

o 

o 

o  o 

o 

o 

O 

O 

O 

o 

O 

# 

. 

in 

ro 

no 

r- 

O 

X  O' 

ro 

X 

X 

<r 

ro 

O  cE  rO 

X 

O  -TO 

n- 

as 

O' 

as 

ro 

ro 

o  rv 

X 

ro 

O  no 

o 

X 

o 

o 

ro 

n- 

O' 

o 

o 

ro  O' 

o 

ro  r- 

X 

o 

r- 

r- 

r- 

r> 

© 

O' 

O 

— 4 

no 

ro 

X 

X 

X 

O 

ro 

ro 

ro  ro 

ro  x  x 

-r 

ro 

ro 

ro 

—4 

O  o-  cd  r- 

ro 

ro 

ro 

O  cc 

X 

o 

no 

n- 

m 

o 

n- 

X 

ro 

rv 

no 

c  © 

© 

cc 

r-~ 

n- 

o 

— 

-• 

— 

no 

ro 

ro 

ro  ro 

ro 

no 

ro 

ro 

ro 

no 

nj 

no 

ro  ro  ro 

no 

ro 

no 

no 

no 

no 

C 

O  O' 

CD 

r- 

r~ 

r~ 

ro 

ro 

ro 

ro 

ro 

ro 

X 

•o 

X 

X 

X 

- 

- 

- 

- 

- 

- 

~ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

~ 

- 

- 

—  - 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

o 

O 

o 

o 

o 

o 

© 

©  o 

o 

o 

o 

O 

O 

o 

© 

© 

ONh 

nj 

o  o 

ro 

r- 

O  ro 

00 

m  ro 

ro 

ro 

ro 

X 

<M  O  ffl 

r~ 

ro  n- 

O  ro  ro 

r-  O' 

as 

o  ro 

n- 

o  ro 

O' 

no 

O  X 

ro 

rv  co 

© 

O' 

no 

ro 

© 

ro  © 

© 

r» 

o 

o ■ 

O 

o  o 

-- 

ro 

ro 

X 

n- 

ro 

ro 

O' 

O' 

O 

o  o  o  o 

O  O  O' 

On  O 

CD 

<D  r- 

ro 

X 

ro 

rv 

rv 

*-* 

O' 

CD 

X 

4-4 

X 

ao 

no 

as 

X 

ro 

O  CD 

© 

r* 

ro 

X 

X 

ro 

ro 

ro 

o 

ro 

<NJ 

rvj 

rvl 

no 

no 

ro 

no 

OJ 

<M 

ro 

ro 

no  ro 

ro  m 

•O 

ro 

ro  rv  ro 

ro 

no 

nj 

ro 

ro 

no 

no 

no  ro 

no 

no 

*“* 

O' 

CO 

CD 

r- 

r~ 

n- 

ro 

ro 

ro 

o  ro 

ro 

ro 

ro 

O 

- 

- 

- 

- 

—4 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

44 

- 

- 

--- 

4-4 

- 

44 

- 

- 

- 

- 

- 

- 

— 

- 

- 

— 

44 

- 

0 

0 

0 

0 

© 

o  © 

o  o 

o 

o 

o 

o  o  o  o 

O  fO 

«T 

O' 

•o 

m 

r- 

X 

o 

ro 

r- 

no 

o  «o 

CD 

ro 

>r 

ro 

®  O  O' 

o 

cd  ro 

ro 

O 

X 

O' 

X 

f" 

o 

r- 

X 

as 

rv 

cr 

O' 

ro 

r> 

ro 

O' 

o  © 

O' 

_ 

o 

no 

© 

o 

x 

x 

in 

m 

■o 

n- 

CO 

ro 

ro 

X 

X 

ro  O 

n- 

r-~ 

n-  n-  n- 

CD 

n- 

r~ 

ro 

ro 

X  X 

ro 

rv 

ro 

O'  as 

ro 

ro 

CD 

nj 

© 

X 

rv 

O' 

© 

© 

r- 

ro  x 

X 

o 

ro 

fo 

ro 

no 

ro 

ro 

no 

ro 

ro 

ro 

fO 

ro 

ro  nv 

ro 

rv 

rv 

ro  ro  ro  rO  ro 

ro 

n 

X 

ro 

X 

m 

ro 

n 

ro 

m  rv 

ro 

•v/w(voff* 

O 

CO  CO 

© 

© 

© 

© 

r- 

r- 

r- 

© 

n- 

r~ 

r- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- - - 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

O 

O 

o 

O 

o 

O 

o 

o 

o 

o 

o 

0 

0 

0 

0 

nj 

© 

rvj 

O  O'  o 

m 

a> 

ro 

o  r~ 

ro 

O  CD 

ro 

O  no  x  CD 

no 

rv 

rv 

ro 

X 

X 

ro 

no  o  n- 

ro 

no 

ro  ro 

r- 

X 

CD 

r- 

ro 

ro  x 

© 

X 

ro 

X 

ro 

ro 

r-  © 

ro 

ro 

o 

c 

O' 

O' 

O 

r*' 

<• 

m 

o 

O' 

O' 

O 

4-4 

—4 

ro 

no 

rv 

rv  rv  rv 

-y 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

rv 

no 

rv 

no 

44 

r- 

n- 

CD 

ro 

O' 

t' 

X 

rv 

ro 

44 

O  © 

r- 

X 

ro 

ro 

ro 

X 

00 

•VI 

no  ro  ro 

r. 

ro 

ro  ro 

ro 

ro 

ro 

ro 

>r 

«T 

«r 

>r 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

ro 

rv 

44 

o 

O  O' 

O' 

O' 

O' 

O' 

o 

O' 

CD 

© 

© 

CD  CD 

© 

cz 

4— 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

m 

ec 

** 

"* 

-* 

•-* 

-4 

44 

44  44  44 

—4 

r~4 

—4 

44 

44 

4—4 

o  ©  o 

o 

o  o 

o 

o 

o 

0 

0 

o 

O 

r- 

m 

•Vi 

OJ 

ff 

>c 

X 

C 

C 

- 9 ’ 

C 

no 

C 

ar 

X. 

n: 

C 

o 

ro 

no 

oar- 

z 

o 

ro 

ro 

XI 

XI 

•o 

n- 

O' 

o 

—4 

no 

•n 

X 

-o 

r- 

ro 

O' 

o 

O 

-•  -•  rsl 

c 

rv 

ro 

ro 

ro  ro 

ro 

ro 

nT 

c 

>T 

>y 

>T 

-T 

•j- 

>3 

X 

XXX 

i/v 

© 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

4-4 

- 

- 

- 

- 

- 

— 

a 

ro 

no 

O 

o 

o 

o 

o 

rv. 

o 

o 

ro 

O' 

nj 

c 

©  rv  rv 

~ 

IT 

cz 

© 

c 

— * 

rv 

r ■ 

m 

r- 

a 

O 

•N. 

ro 

X 

ro 

ro 

cr 

c-  c  — 

>-  r- 

ro 

ro 

ro 

ro 

ro 

ro 

'T 

'T 

xv 

X 

X 

X 

X 

X 

X 

X 

x  x  ro  ro 

QC 

- 

4— 

- 

- 

44 

- 

- 

- 

- 

- 

4-4 

- 

- 

4-4 

4-4 

—4 

- 

- 

4-4  4^-* 

*e 

no 

© 

r 

ro 

Xl 

CO 

fs*  cc 

ro 

ro 

nj 

no 

>T 

ro 

rv 

a. 

—  roc 

o  ro 

uv  in  O  N 

CD  O' 

—• 

ro 

o  r-  ro  O' 

O 

44 

ro 

-r 

x  n-  co  o 

ro 

ro 

ro 

ro 

>T 

ro 

-T 

X 

in 

XV 

X 

X 

X 

X 

ro 

ro 

ro 

ro 

ro 

ro  ro  r- 

- 

44 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

— 

- 

- 

- 

- 

- 

--- 

— *  in 

rO 

O' 

ro 

•r 

ro 

*— 4 

_ 

ro 

■•j 

44 

r- 

x  n- 

ro  no  o 

c  c 

© 

O' 

©  — •  ro 

>c 

cr 

c 

ro 

fO 

X 

n- 

a 

c 

»4 

X 

ro  cd 

ro 

ro 

ro 

ro 

m 

m 

m 

m 

m 

in 

-o 

O 

>o 

vO 

o 

ro 

ro 

n- 

r- 

r~  r~ 

- 

- 

- 

- 

- 

- 

—* 

- 

- 

- 

4-4 

- 

- 

- 

- 

- 

->4 

- 

--- 

•  IT 

j 

x 

r- 

© 

v> 

•c 

c 

n. 

05 

r- 

ro 

n~ 

r 

X 

—  c  e 

o  rv 

m 

o 

x  O' 

ro 

in 

O 

O  ro 

X 

r- 

cn 

O  ro  rv  X 

m 

IT 

x 

IT 

X 

m 

m 

X 

«C 

o  <c 

-c- 

-c 

r~ 

r- 

r- 

n> 

n- 

CD 

ec  ec  x 

- 

- 

44 

- 

*-• 

•—4 

- 

- 

- 

- 

- 

- 

- 

— 

- 

- 

- 

4-4 

ro 

o 

® 

O  ro 

X 

ro  o  o  in 

o 

o 

r- 

o  ro  no  o  n- 

C 

c 

4— 

rv 

X 

c: 

O 

rv 

X 

c 

. — 

no 

-c 

ro 

C 

no 

•c 

X 

r*  o  c 

Csl  O 

o 

•O 

■o 

o 

o  n- 

r- 

r- 

r- 

ro 

ro 

ro 

ro 

ro 

o 

o 

c> 

O'  o  o 

~4  44  ro 

"* 

rv 

r 

n- 

a 

O' 

o 

O' 

O  X 

_ 

c 

X 

O' 

ro 

ro 

n- 

a  n- 

ro  ro  O' 

O  f- 

X 

rvi 

X 

o 

■r 

r- 

3 

X 

n- 

O' 

ro 

n-  O' 

O  no  st  ro 

— 

<C 

< 

-c. 

r- 

n- 

r- 

X 

c 

a 

C 

c 

c 

c 

o 

C 

C 

c 

c 

—  4-4  4-4 

J, 

J 

Ji 

J 

J, 

— 

J, 

^4 

Ji 

^4 

J. 

— 

— 

4- 

ro 

ng 

no 

ro 

no 

ro  no  ro 

X  o 

o 

o  o  o 

o  o 

o 

o  © 

©  o 

© 

o 

o 

o 

©  x 

©  x  o  x 

O  x 

o  x  o  x 

oxoxoxoxox 

o  x 

o  x 

©  X 

O  X 

©  x 

o 

©  o 

o 

o 

o  ro 

ro 

ro  O  © 

ro  ro 

ro 

o  © 

ro  x 

n* 

rv 

no 

o  © 

O  CD  CL  O- 

r-  ro 

«ro  x  x  ro 

ro  rv  m  nn  -  -o 

©  © 

cd  r- 

©  ro 

ro  x 

x  ro 

ro 

rv  ro 

© 

r  ro 

ro 

ro  ro  r 

rv  r 

r 

rv  oj 

rv  no 

no 

ro 

ro 

no 

no  44 

44-4-4- 

4—  44 

—  4-  44  — 

c  c 

c  c 

c  © 

c  c 

C  C 

c 

c  c 

c 

© 

Elil 8  '.v,;.'.  i 
Jfi\  jtt.Our  J__-* 


6 


XO  x/i)0 


Figure  3.  1977  ROBIN  Program:  Drag  Table 


OhAG  TAhll  11AUG  1977) 


) 


I 


L 


c 


IT 

s~ 


z 3  x  c* 

3  er  in.  y  y  p  fM  im  rv  3  >  O'  x 


•  -  n  >c  m  «  r  :  £  n-it  c 

C  -  3  3  X  P 

C  •nyyyyyyyyy-''-''  — 


n  cccocccccccc-c 


in*-^cxC'PfNy©er'*,<N©  c 

T'  "■  o  c  ir  #  c  i*'  \  -  3  c  x 

mininyyyyyyyypppp. 


C333C33CC  3  C  C  s  c  3 


• 

— * 

S3 

«o 

>r 

P 

o 

>r  — 

>T 

P 

O' 

O 

©  -y 

P  © 

4 

3 

C 

<. 

X 

— 

o 

©  n\ 

p 

IN 

IN. 

— 

j' 

X 

p 

p 

o 

X 

P 

© 

n 

m 

m 

>r 

>r  -r 

>r 

vT 

'T 

>r 

>r 

p 

p 

p 

p 

y 

3 

3 

3 

3 

C 

C 

3  © 

c 

3 

© 

e 

c 

3 

C 

3 

C 

. 

in 

© 

© 

rsi 

O' 

c 

^  -T  3 

X. 

— 

r- 

IT 

c 

_ 

in.  O' 

r- 

c 

—  fpl 

fTN 

O' 

C  p 

X 

X 

- 

_, 

_ 

*n 

IN 

O 

—• 

— • 

p 

c 

<4 

C  O  C 

O' 

O' 

<n 

X 

P- 

X) 

•o 

m 

>P  »NI  — 

O' 

p- 

in 

>r  — 

X 

in  o 

m 

— 

O' 

p  m 

>r 

p 

P 

IM 

—4 

rp 

® 

*4 

P 

S\ 

C 

c 

c 

C 

c 

c 

C  C  C 

O' 

c 

O' 

O' 

e 

o 

o 

O' 

O'  c  o 

X 

X 

c 

x  or 

p 

*- 

m 

ir. 

© 

<r 

■n  >r 

N* 

•n 

*n 

' T 

r- 

P 

rr 

- 

- 

- 

- 

- 

- 

-  -  - 

« 

c 

= 

c 

= 

= 

- 

C 

o  c  s 

- 

c 

- 

3  3 

3 

- 

3 

= 

= 

=■ 

3  C 

=- 

c 

= 

•i 

= 

- 

m 

© 

m 

fM  O' 

o  «r  o 

•c 

>n 

cc 

INI 

f*\  O  O 

3 

INI  P- 

m  x 

© 

p 

INI 

X 

in 

O' 

P  IM 

c 

p 

•T 

in 

© 

fM 

- 

3 

fM 

IN 

TNI 

rvj 

04 

'VJ 

4)  O  M 

o 

coo 

c 

r- 

■C  «r\  >r  ini 

— 

O' 

O  o  INI 

O' 

o 

3 

in 

IM  o 

P  © 

n 

>r 

P 

IM 

— 

3 

O' 

<r 

X 

c 

O 

o 

© 

O 

c 

c 

CSC 

o 

O  O' 

O' 

O'  O' 

o 

O' 

O'  o  o 

O' 

X 

a  x 

P 

© 

© 

•n 

m 

>r 

>r  >r 

'T 

>r 

'T 

y 

y 

p 

P 

P 

p 

- 

- 

- 

- 

- 

- 

-  — 

- 

- 

c 

o 

© 

© 

O 

c 

03  0 

o 

o 

C 

3  O 

o 

<= 

C 

® 

o 

o  © 

C 

c 

c 

3  ©  C 

o 

3 

o 

. 

X 

O' 

p 

•T 

7\J 

rj 

o*  m  m 

o 

>r 

S'. 

S  <N 

O  x  in 

O' 

c 

X 

c  o 

- 

c 

O'  o 

© 

in 

© 

c 

p 

r 

3 

p 

p 

p 

m 

m 

7»> 

<M  rw  rg 

—• 

*-4 

— 

O  O  O' 

x  h- 

o  «n  -r 

i\l 

— 

CD 

p* 

O  p 

— * 

© 

p 

3  P  © 

m 

>r 

p 

INI 

—• 

3  O'  O' 

X 

C 

c 

© 

C 

C 

c 

o 

c  c  c 

c 

c 

c 

c 

C  O' 

O' 

O' 

O'  O'  O' 

O' 

O' 

X 

X  X 

® 

© 

© 

m 

m 

m 

>r  >r 

vT 

'T 

y 

y 

p 

P 

P 

im 

- 

- 

- 

- 

- 

- 

--- 

- 

- 

- 

- 

- 

c 

C 

c 

o  o  c 

C 

= 

c 

©  o 

c 

= 

© 

© 

c 

= 

3  3 

3 

o 

o 

C 

c 

3 

c 

C 

3 

p 

© 

sT 

r«j  7^7  O  X)  oj 

X 

O'. 

O  r- 

sT 

f\i  «r 

r-  r-  o 

O  ^ 

o 

—  -r 

® 

© 

C  IT- 

m 

x. 

p  ^ 

C 

c 

fM 

in 

© 

y 

in 

P 

P 

3 

y 

v» 

-r 

w  O 

<NJ 

INI 

INI 

—4 

— • 

3  C 

X 

p-  -c  n 

>T 

IM 

3 

x  «n 

•-4 

X 

iM  P  P 

3 

X  •" 

© 

in 

-S' 

p 

INJ 

3 

c 

X 

c 

c 

c 

o  e 

o 

3 

o  c  c 

3 

c 

c 

o 

c 

O  O' 

O' 

O'  O'  O' 

O' 

O' 

o* 

®  X 

X 

© 

© 

n 

S' 

-t  >r 

•n 

>r 

*n 

y 

y 

y 

p 

P 

<M 

- 

- 

- 

- 

- 

- 

-  -  - 

- 

- 

- 

- 

- 

- 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

0 

0 

0 

0 

© 

o  o  o 

o  o 

0 

0. 

0. 

o  o 

o 

, 

P 

-vT 

fM 

c  x 

X  rvi  O' 

o 

r- 

f\l 

X 

C  fp 

•£ 

c  r-  x 

®  O' 

o 

im  <r 

O' 

© 

O' 

r-. 

®  p 

p 

IM 

p 

O' 

c 

X 

O' 

o 

m 

i© 

© 

»r* 

>r 

-r  m 

ro 

INI 

INI 

—< 

•-4 

O  O 

co  -o 

-T 

IN) 

c 

O'  o 

INI 

O' 

IM  P 

>r 

^4 

CD  P 

vT 

p 

fM 

fM 

O' 

® 

9 

O 

C 

<3 

o 

o 

z 

C  S  O 

o 

o 

a 

3 

o 

o 

O  O 

O'  O'  O' 

O' 

O' 

c 

x  ® 

X 

© 

© 

n  s\ 

n 

*r  -r 

>r 

>T 

•y 

y 

y 

y 

p 

p. 

- 

- 

- 

- 

- 

- 

--- 

- 

- 

- 

- 

- 

- 

3 

0 

0 

0 

C 

c 

3 

C  3 

3 

3 

c 

3 

c 

© 

©  O 

c 

c 

3 

c 

3 

3 

c 

c 

c 

. 

P 

y 

p 

rr\ 

r^v 

-*C  -f 

r- 

•Q 

_ 

>r 

S3 

•o 

—  — •  O 

c 

INI  O' 

P»  i\i 

m 

® 

O' 

® 

c  <■* 

O  O' 

ini  in 

P 

® 

O' 

© 

© 

c 

© 

© 

X) 

<o 

O 

•C  *rs  in 

S' 

m 

••A 

r-» 

INi 

— 

—4 

O  O'  ® 

o 

—• 

C  X 

>T 

c 

® 

m 

*-4 

O'  X 

p 

n 

y 

p 

IM 

—< 

C  O' 

c 

c 

c 

c 

O 

c 

C  C  C 

c 

o 

C 

c 

o 

o 

c 

C  O'  O' 

O' 

O' 

O' 

O'  CO 

X 

p 

© 

in 

m 

»n 

©  'T 

nT 

V* 

•y 

«y 

y 

y 

y 

y 

p 

-* 

- 

- 

- 

- 

- 

- 

-  — 

- 

- 

- 

- 

- 

- 

— 

— 

-oo 

o 

c 

© 

o  © 

o 

c 

o 

o 

3 

C 

©  o 

o 

o 

c 

c 

=> 

3 

© 

c 

c 

. 

25 

_ „ 

P 

m 

_ 

S' 

IM 

© 

IN 

*M 

o  « 

-r  ^  r> 

■G 

•c 

f'J 

p* 

C 

—  CO 

«_ 

•J- 

X  w 

p 

m 

INi 

r 

iM 

c 

©  © 

O' 

p 

X 

n 

© 

3 

-5 

-2 

y 

v?’ 

<n 

m 

© 

r- 

X 

X 

X 

X 

r-  r- 

o 

o 

n 

•T 

vT 

Ni 

INI 

—  C  X 

p- 

m 

er\ 

—  O' 

n 

—4 

m 

3 

© 

IM 

©  ® 

p- 

© 

n 

y 

p 

IM 

~ 

~ 

c 

3 

- 

3 

z 

C 

3 

- 

~ 

3 

© 

c 

c 

- 

c 

3  c  c 

s 

- 

- 

3 

c 

- 

3  3  C* 

O' 

O' 

O 

O'  C3 

X 

P 

© 

© 

m 

in 

n  vf 

V* 

s 

©■ 

y 

y 

y 

y 

y 

- 

- 

- 

- 

— 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

--- 

- 

- 

- 

— 

— 

- 

— 

—  —  C 

o 

3 

o 

c  © 

© 

3  © 

o  c 

O 

O  O 

© 

© 

o  o  c 

c 

C 

© 

© 

m 

<n 

in  in 

X 

IM  © 

O' 

(M 

y  © 

p 

© 

*n 

t*s 

r-  <b  -o 

O' 

<c 

vT 

O 

>y 

P- 

P-  -O  INJ 

r- 

<c  c 

O' 

P 

P 

m  ® 

O' 

c 

IM 

«n 

p 

in 

m 

•n 

y 

3 

n 

m 

j" ■ 

n 

lT 

© 

SC. 

— 

p 

P 

p 

O' 

O' 

j 

O' 

co  x  x 

e 

c 

•C 

S' 

V* 

-T 

i*> 

m  -*  C 

X 

c. 

>r 

r*\  — 

c 

IN  © 

—4 

p 

r- 

—  3 

® 

p 

© 

in 

p 

IM 

p 

- 

c 

c 

©  o 

c 

- 

© 

o 

c 

o 

c 

c 

c 

O 

c 

c 

o  c  c 

C 

c 

D 

Z> 

3 

c* 

c 

O  C  O  O' 

O' 

O' 

O'  o 

X 

O 

© 

in 

m 

*n  -n 

>r 

>r 

>r 

y 

y 

y 

y 

y 

- 

- 

- 

— * 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

-  -  - 

- 

— 

— 

— 

- 

— 

~ 

— 

—  -  — 

o 

o  © 

c  © 

c 

©  O 

o  © 

O 

coco 

o 

0 

‘0 

‘0 

O  © 

© 

# 

> 

i> 

p 

in 

r 

n 

y 

i*-* 

r> 

n 

x 

o 

, 

C  O' 

X 

m  t*\  t»j 

in 

r*> 

in 

^  J  A. 

*n 

o 

p-  -n 

. 

c  ® 

m 

p 

p>  n 

3  P 

n 

p 

y 

y 

p 

n 

c 

■o 

© 

©  o 

p 

r- 

CO 

P 

p 

04 

o4 

O  c 

C  C  C 

c 

O' 

X 

p*- 

P- 

X) 

o  m 

P\  IN, 

3  « 

o 

>r  ini 

O' 

n 

p 

INI  X 

>r 

INI  — 

C  3 

® 

p 

©  »n 

y 

p 

p 

p 

3 

r 

c 

3 

C 

3 

C 

~ 

— 

— 

— 

— 

— 

— 

—  —  — 

c 

c 

' 

c 

- 

C 

c 

c  c  c 

o 

O' 

O'  o- 

X 

p- 

© 

© 

m 

m 

m  m 

in 

-J- 

y 

y 

y 

y 

y 

— 

— 

— 

— 

— 

— 

- 

— 

— 

- 

- 

- 

— 

- 

p 

— 

- 

- 

—  —  - 

— 

- 

- 

- 

- 

- 

— 

— 

—  —  — 

- 

C 

3 

C  3 

c 

3 

c 

3 

C 

c 

3  3 

3 

3 

3 

C 

3 

3 

c 

3 

3 

— 

© 

© 

©  o 

3 

3 

O 

o 

o 

c 

3  © 

© 

o 

3 

C 

o  in 

©  n 

C  in  ©  m 

On  on 

o  in 

O  in 

©  in  ©  n 

©non 

©  in 

O  n 

©  in 

O  n  o 

© 

©  o 

o 

» -> 

© 

y 

IM  O 

T) 

© 

y 

IM 

c 

X 

©  m 

y 

p 

<M 

*“• 

O  O' 

O'  ® 

®  p  p  © 

©  n  in  y 

y  p 

P  <M 

INI  p  —  o 

©  O'  O'  ® 

®  p 

P  © 

©  m 

n  y  y 

p 

fM  — 

o 

T 

y 

y 

y  y 

P 

P 

p 

P 

p 

fM 

<M  IM 

IM 

IM 

IM 

fM 

fM  — • 

—  — 

p  p  p  - 

-  -  -  - 

p  - 

p  p 

p  p  p  -1 

—  03© 

©  o 

C  O 

3  © 

o  ©  c 

3 

c  © 

c 

7 


Figure  3  (cont'd).  1977  ROBIN  Program:  Drag  Table 
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Figure  3  (cont’d).  1977  ROBIN  Program:  Drag  Table 


Comparison  of  Output 


Comparison  of  a  flight  reduced  with  both  the  1977  ROBIN 
and  the  1972  HIROBIN  program  is  shown  in  Figure  4.  Note  the 
difference  in  the  initial  values  of  temperature  at  97  Km 
(due  to  the  1976  Standard  Atmosphere) ,  the  gradual  change  in 
temperature  values  above  76  Km  (removal  of  supersonic  bias) , 
the  removal  of  the  Mach  1  spike  near  68  Km,  and  the  drag  table 
induced  temperature  changes  below  58  Km. 
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Fiqure  4.  Comparison:  1977  ROBIN  vs.  1972  IIIROBIN 


SECTION  2 

PURPOSE  AND  USE  OF  PROGRAM 


The  1977  ROBIN  program  is  specifically  designed  for  computer 
reduction  of  ROBIN  flights  which  attain  an  apogee  of  approximately 
115  Km,  and  are  tracked  by  an  FPS/16  radar.  The  program  also 
provides  accurate  reduction  for  ROBIN  flights  attaining  other 
apogees  between  65  and  150  Km.  If  a  passive  sphere  other  than 
the  one  meter  ROBIN  (approximately  165  grams)  is  used,  reduction 
with  this  program  is  still  possible.  There  may  be  times,  however, 
when  drag  coefficients  are  not  available  and  thus  density, 
temperature,  and  pressure  cannot  be  computed.  The  time  of  fall 
test  used  to  determine  sphere  collapse  is  not  used  for  spheres 
with  dimensions  that  differ  from  the  ROBIN. 

At  each  second  of  flight,  the  program  calculates:  altitude, 
wind  components,  and  direction;  density,  temperature,  and  pressure 
the  RMS  noise  error  in  winds,  density,  temperature,  and  pressure, 
and  other  balloon  dynamics. 

The  following  are  some  of  the  important  features  of  the 
1977  ROBIN  program. 

(1)  The  program  begins  wind  and  thermodynamic 
computations  when  the  vertical  acceleration  (*z*) 
becomes  greater  than  -8  m/sec2.  For  a  sphere 
apogee  of  115  Km,  this  occurs  at  about  95  Km. 

(2)  The  program  computes  the  RMS  noise  errors 
in  wind,  density,  temperature,  and  pressure.  A 
discussion  of  other  type  errors  is  contained 

in  Section  4  of  this  manual. 

(3)  The  program  determines  sphere  apogee  whenever 
sufficient  tracking  data  is  available. 

(4)  The  program  uses  two  tests  as  an  indication 
of  sphere  collapse,  a  time  of  fall  test,  and  a 
density  gradient  test  (X  check) . 

(5)  At  the  end  of  the  flight  the  program  prints  a 
frequency  response  table  for  density  and  winds. 

This  enables  the  user  to  observe  what  wavelengths 
have  been  retained  in  the  output  data. 
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(6)  The  program  provides  the  following  options 
of  output:  a)  output  at  each  second  of  flight; 

b)  additional  interpolated  output  at  each  200  meters 
of  altitude;  and  c)  output  each  kilometer  only. 

(7)  Below  60  Km  the  density  smoothing  interval 
expands  in  such  a  manner  so  as  to  maintain  a 
constant  two  percent  noise  error  in  density. 

(3)  If  spheres  other  than  the  ROBIN  are  used, 
the  reduction  can  be  optimized  by  changing  the 
smoothing  input  parameters.  However,  when  using 
a  ROBIN  sphere  at  any  escape  altitude,  it  is  of 
fundamental  importance  that  one  not  change  any 
of  the  smoothing  input. 
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SECTION  3 

OPERATING  INSTRUCTIONS 

The  High  Altitude  ROBIN  program  was  written  using  FORTRAN  IV. 

The  program  was  tested  on  the  CDC  6600  computer  at  Wright-Patterson 
Air  Force  Base. 

The  basic  program  generally  requires  two  to  three  minutes 
to  run  on  the  CDC  system;  however,  the  run  time  can  be  almost 
twice  that  if  the  data  tape  requires  a  great  deal  of  editing. 

The  input  to  the  ROBIN  program  consists  of  time,  X,  Y  and  Z 
coordinates.  Time  is  in  units  of  seconds  Zulu  and  the  X,  Y  and  Z 
position  coordinates  in  units  of  feet.  Generally  the  position 
coordinates  are  given  with  respect  to  the  coordinate  system  of  the 
firing  range.  The  output  from  the  ROBIN  program  gives  wind 
measurements  with  respect  to  the  north  and  east  components.  To 
accomplish  the  transformation  of  coordinate  systems  the  positive 
clockwise  angle  from  geographical  north  to  the  positive  X-axis  of 
the  range  coordinate  system  must  be  known.  The  angle  is  called  Z^ 
and  is  input  to  the  ROBIN  program. 

The  ROBIN  system  is  divided  into  two  related  programs. 

The  first  is  the  tape  editing  program.  It  should  be  ncted  that 
this  program  has  not  changed  since  the  1972  version.  The  main 
function  of  this  program  is  to  replace  bad  data  points  with 
interpolated  data.  Its  printed  output  consists  of  the  discarded 
old  points  and  the  added  new  points.  (See  Figure  5).  The 
editing  routine  then  writes  the  edited  data  on  a  tape  or  disk. 

The  data  is  written  in  97  word  records  with  the  exception  of  the  first 
record  on  a  file  which  contains  the  flight  identification  number. 

The  97  word  records  begin  with  a  FORTRAN  loader  word  and  thereafter 
contain  a  96  word  array  of  data  points.  An  undetermined  number  of 
records  are  written  per  file  since  a  new  file  is  started  every 
time  a  gap  greater  than  two  seconds  occurs  in  the  data.  Each  time 
an  end  of  file  is  written  a  notice  is  placed  in  the  printed  output. 


DATA  INPUT  TO  ROBIN  EDITING  PROGRAM 


Card  1 

Column 

Variable 

Description 

1-10 

ID 

Flight  Identification  number 

11-15 

IX 

Output  unit 

15-20 

IY 

Input  Unit 

21-30 

XTOL 

The  X  tolerance 

=  150.0  ft. 

31-40 

YTOL 

The  Y  tolerance 

=  150.0  ft. 

41-50 

ZTOL1 

The  Z  tolerance 
ft.  =  360.0  ft. 

at  or  above  200000 

51-60 

ZTOL2 

The  Z  tolerance 
=  200.0  ft. 

below  200000.0  ft 

61-70 

ALT 

Read  past  ascent  data  to  ALT  = 

300000.0  ft.  (for  high  apogee  Robin; 
adjust  if  apogee  is  known  to  be 
lower ) 


The  main  FOB IN  program  uses  the  edited  data  to  compute  the  de¬ 
sired  meteorological  parameters.  The  methods  used  to  compute  the 
parameters  are  discussed  in  References  2,  4,  and  5.  Once  the  parameters 
are  computed  they  are  printed.  The  ROBIN  program  provides  the 
options  of  printing  only  the  actual  data,  printing  the  actual  data 
and  data  interpolated  at  every  200  meters  of  descent,  or  printing 
interpolated  data  every  kilometer.  Figures  6,  7,  and  9  illustrate 
these  alternatives. 

DATA  INPUT  TO  ROBIN  PROGRAM 

Card  1 


Column 

Variable 

Description 

1-5 

NXY1 

Number  of  points  for  wind  velocity 
fit  *  51. 

6-10 

NXY2 

Number  of  points  for  wind  acceler¬ 
ation  fit  *  35. 

11-15 

NZ1 

Number  of  points  for  density 
velocity  fit  =  19. 

16-20 

NZ2 

Number  of  points  for  density 
acceleration  fit  =•  21 

21-25 

IDEGXl 

Degreee  of  fit  for  wind  velocity=3 

26-30 

IDEGX2 

Degree  of  fit  for  wind  acceleration 
-  3. 

31-35 

IDEGZ1 

Degree  of  fit  for  density  velocity 
-  1. 

36-40 

IDEGZ2 

Degree  of  fit  for  density  acceler¬ 
ation  *  3. 

Card  2 

1-80 
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Description  title  of  job. 

Card  3 
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ID 

Flight  number. 
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IX 

Input  tape  (disk)  unit. 
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IY 

Output  tape  (disk)  unit. 
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Latitude  of  launch  site 
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Height  above  mean  sea  level  of 
launch  site  (meters) . 

28-33 

ZB 

Angle  of  X  axis  from  north  (degrees) . 
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Figure  6.  University  of  Dayton  1977  ROBIN  Program 
(Actual  Data) 
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Figure  7.  University  of  Dayton  1977  ROBIN  Program  (Data 
Interpolated  every  200  meters) 
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/ 


X 


34-40 

41-48 

49-51 


AMS 


Mass  of  Balloon  (kilograms) . 

DIA  Diameter  of  balloon  (meters) . 

Interp  Switch  Interp  ^  0  no  interpolation 

Interp  =>=  1  interpolation  every 
200  meters.  Interp  =  2  print  at 
every  kilometer  only. 

The  deck  setup  is  shown  in  Figure  9. 


Input  and  Output  Units  by  ROBIN  System 
Units  five  and  six  are  the  standard  FORTRAN  input  and  output 


units . 

Units  Mode 

Edit  Program 

Block  Size 

Usage 

8 

Binary 

97 

Edited  input  tape  (or  disk 
output  for  edit  program. 

1 

Binary 

97 

Input  tape 

ROBIN 

Program 

8 

Binary 

97 

Edit  input  tape  (or  disk) . 
Same  as  unit  8  in  edit 

program. 

10 

Binary 

6 

Output  tape (or  disk)  for 
plot  program. 

20 


Figure  9.  Typical  Deck  Setup 
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SECTION  4 

DATA  INTERPRETATION 

Apogee  Determination 

Sphere  apogee  can  be  determined  if  the  radar  acquires  track 
before  the  downward  velocity  of  the  sphere  exceeds  150  m/sec. 

The  sphere  apogee  must  be  determined  in  order  to  use  the  time  of 
fall  test  between  80  and  55  Km.  If  radar  acquisition  does  not 
permit  apogee  determination,  the  printout  "Balloon  Apogee  Unknown" 
appears  and  the  time  of  fall  test  of  sphere  collapse  commences 
at  55  Km. 

Sphere  Collapse 

a)  Time  of  Fall  Test 

The  estimated  time  of  fall  of  a  ROBIN  sphere  through  5  and 
10  Km  altitude  layers  in  the  1962  Standard  Atmosphere  density 
(p„62„)  and  +  20  percent  variations  in  density  from  the  Standard 
Atmosphere  is  given  in  Table  1.  If  the  actual  time  of  fall  of  a 
sphere  is  not  between  that  estimated  for  a  0 . 8 p „ 62  „and  1.2pll62„ 
density  profile,  the  message  "Time  of  Fall  Test  Between  XX  and 
XX  Km  indicates  Collapse"  is  printed  out.  Invariably  when  a 
sphere  collapse  occurs,  the  actual  time  of  fall  greatly  exceeds 
the  imposed  tolerances.  If  the  time  of  fall  test  is  satisfied,  the 
message  "Balloon  Still  Inflated  at  XX  Km"  appears  in  the  printout. 

b)  Density  Gradient  Test  (X  Check) 

Below  60  Km  the  density ■  gradient  (  — lean  be  estimated 

\  P  dz  / 

by  computing  A  defined  as: 

X  -  (Lo‘3(i2>  -  ‘•“'J'v)  *  l  5§ 

The  expected  value  of  A  between  30  and  60  Km,  obtained 
from  the  "62"  Standard  Atmosphere,  is  A  =  .00014.  To  allow  for 
density  perturbations  in  the  real  atmosphere,  a  tolerance  of  .00005  is 
allowed.  In  addition  to  allowing  a  tolerance  in  A  for  real  density 
perturbation,  a  tolerance  must  also  be  allowed  for  the  inaccuracies 
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Apogee  (km) 


100  km 


80  km 


70  km 


p62 

1.2p,, 


Table  1.  Time  of  Fall  in  Seconds  through  Standard 
Atmosphere  and  +  20%  Deviations  from 
Standard  for  Sphere  Apogees  Between  70 

ana  .50  Kin 


of  the  parameters  used  to  compute  A .  For  a  two  percent  noise  error 
in  density,  the  3o  noise  tolerance  for  A  is  3o^  =  • 

The  check  is  summarized  as  follows: 

Expected  value  of  A  =  .00014; 

Tolerance  for  variation  in  density  is  +  .00005;  and 

08484 

3o  noise  tolerance  is  +  - — -  ,  Az  is  in  units  of  meters. 

-  Az 


A  sphere  is  considered  inflated  if: 


00014 


-[<• 


00005)  +( 


.08484 

Az 


] 


1/2 


2  08484  ^ 

^  A  <  .00014+  (.00005)  +(— 


If  the  A  check  fails,  the  message  "Density  Gradient  Exceeds 
Nominal  Value"  appears  in  the  printout.  The  A  check  is  the  secondary 
balloon  collapse  check  and  should  be  used  in  conjunction  with 
the  time  of  fall  tests.  The  A  check  is  a  sensitive  check  and 
occasionally  can  be  tripped  by  a  vertical  wind  or  density  anomaly 
not  connected  with  sphere  collapse.  In  this  case  the  error  message 
will  appear  only  sparingly  and  the  time  of  fall  test  will  validate 
proper  balloon  inflation.  When  sphere  collapse  is  indicated  by 
the  time  of  fall  test  over  say  a  5  Km  interval,  the  A  check  will 
often  pinpoint  the  actual  collapse  within  an  accuracy  of  1  Km.  If 
a  non-ROBIN  sphere  is  used,  the  A  check  serves  as  the  only  indica¬ 
tion  of  sphere  collapse. 

Density  Errors 

Sources  which  can  produce  significant  density  errors  are 
as  follows. 

a)  Random  error  in  density  due  to  radar  noise. 

b)  Bias  error  in  density  introduced  by  smoothing  function. 

c)  Error  in  initial  guess  at  temperature. 

d)  Drag  coefficient  errors. 

e)  Vertical  winds. 


a)  Noise  Error  in  Density 

The  random  error  in  density  is  that  error  resulting  from 
noise  in  the  radar  signal.  The  magnitude  of  this  error  for  a  nom¬ 
inal  115  Km  apogee  and  the  19-21  linear  cubic  smoothing  (expanding 
to  51-21)  used  on  density  is  given  by  the  first  line  of  Table  2. 

The  two  percent  noise  error  in  density  is  maintained  below  60  Km  by 
gradually  expanding  the  smoothing  interval  from  19-21  to  51-21 
linear  cubic.  The  RMS  noise  error  in  density  is  computed  in  the 
ROBIN  program  by  assuming  an  RMS  error  of  0.15  mil  in  azimuth  and 
elevation  angles  and  6  meters  in  slant  range.  The  RMS  error  value 
printed  out  by  the  program  may  differ  slightly  from  that  given  in 
Table  2  due  to  a  different  sphere  apogee  and  trajectory. 

b)  Bias  Errors  in  Density 

A  bias  error  is  introduced  into  density  when  the  smoothing 
function  does  not  fit  true  perturbations  in  the  position  coor¬ 
dinates.  Incorporated  into  the  1977  program  is  a  bias  removal 
technique  that  will  remove  the  bias  error  resulting  from  a  sphere 
flight  through  the  1976  Standard  Atmosphere.  Since  an  actual 
atmosphere  will  differ  somewhat  from  the  "76"  Standard,  not  all 
of  the  bias  will  be  removed.  The  second  row  of  Table  2  shows 
typical  values  of  the  unremoved  bias  based  upon  simulated  flights 
through  other  atmospheres. 

c)  Initial  Guess  of  T 

- o 

A  guess  of  the  temperature  at  the  initial  altitude  of 
thermodynamic  computations  (approximately  95  Km)  is  required.  This 
temperature  error  produces  a  density  error  in  succeeding  computa¬ 
tion  which  decreases  with  altitude.  The  percent  error  in  density 
produced  by  a  ten  percent  error  in  the  initial  temperature  at  95  Km 
is  given  by  Table  2.  The  effect  of  this  error  becomes  negligible 
after  5  Km  of  flight. 

d)  Drag  Coefficient  Errors 

Errors  in  the  drag  coefficient  effect  density  accuracy  in 
two  ways:  a)  drag  coefficient  error  directly  produces  density  errors 
since  CQ  is  in  the  formula  used  to  compute  density;  and  b)  an 
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Table  2.  Percent  Error  in  Density  from  Various  Sources 
for  Sphere  with  Apogee  115  Km. 


error  in  CD  produces  a  density  error  which  in  turn  produces  an 
error  in  Reynolds  and  Mach  numbers  so  that  the  Reynolds  and  Mach 
numbers  used  to  determine  CD  puts  one  in  the  wrong  region  of  the 
drag  table. 

The  drag  table  used  in  the  1977  ROBIN  program  is  based 
upon  Bailey  and  Hiatt's  (Reference  1)  ballistic  sphere  measure¬ 
ments  with  revision  in  the  subsonic  region  recommended  by  the 
Datasonde-sphere  comparisons  study  of  Engler  and  Luers  (Reference 
2).  Bailey  and  Hiatt  quote  a  maximum  error  of  two  percent  through¬ 
out  the  observational  range.  With  the  revisions  due  to  comparison 
with  other  sensors,  it  is  believed  that  the  drag  table  accuracy 
in  some  regions  is  greater  than  two  percent.  Table  2  gives  the 
best  estimate  of  density  errors  resulting  from  inaccuracies  in 
the  drag  table. 

e )  Vertical  Winds 

In  order  to  compute  density  from  a  passive  sphere,  one 
must  assume  no  vertical  winds.  Actual  vertical  winds  will  produce 
fictitious  density  perturbations.  The  percent  error  in  density 
produced  by  a  1  m/sec  vertical  wind  of  wavelength  2  to  10  Km  is 
given  in  Table  2.  To  obtain  the  percent  density  error  produced 
by  vertical  winds  of  some  other  magnitude,  say  X  m/sec,  the  entries 
in  Table  2  are  multiplied  by  X. 

Pressure  Errors 

Pressure  is  computed  in  the  ROBIN  program  by  an  integration 
of  density.  Hence  errors  in  pressure  can  result  from  all  sources 
of  density  error. 

a)  Noise  Error  in  Pressure 

The  percent  error  in  pressure  resulting  from  the  noise 
errors  in  density  is  given  in  Table  3.  This  error  is  printed  in 
the  computer  output.  It  may  differ  slightly  from  the  values 
given  in  Table  3  due  to  a  different  sphere  apogee  and 
trajectory. 

b)  Bias  Errors  in  Pressure 

The  percent  bias  error  in  pressure  resulting  from  the  un¬ 
removed  bias  in  density  is  given  by  line  2  of  Table  3. 
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The  bias  in  pressure  due  to  a  ten  percent  error  in  the 
initial  temperature  is  given  by  line  3  of  Table  3. 

Error  in  pressure  due  to  drag  coefficient  errors  or 
vertical  winds  cannot  be  directly  calculated.  Since  pressure 
is  computed  by  integrating  density,  the  error  in  pressure  resulting 
from  CD  errors  depends  upon  the  entire  error  profile.  If  CD 
errors  are  biased,  say  all  positive,  the  effect  on  density  error 
will  be  greater  than  if  the  same  magnitude  CD  errors  are  randomly 
positive  and  negative.  Since  the  type  CD  error  profile  is  beyond 
the  present  state  of  knowledge,  pressure  errors  resulting  from 
CD  errors  cannot  be  estimated. 

The  effect  of  vertical  winds  on  pressure  accuracy  depends 
upon  the  vertical  wind  profile  of  the  entire  atmosphere.  Since 
this  information  is  also  beyond  present  knowledge,  the  pressure 
error  produced  by  vertical  winds  cannot  be  estimated. 


Temperature  Errors 

Temperature  is  computed  by  the  gas  law  from  density  and 
pressure.  Hence,  the  sources  of  temperature  errors  are  all  the 
sources  of  density  errors.  The  relationship  between  temperature, 
density,  and  pressure  errors  is  indicated  by  the  formula: 


dT  _  dp  _  dp 
T  "  p  p 


(1) 


Hence,  the  temperature  error  profile  is  obtained  from  the 
density  error  profile  and  the  pressure  error  profile  produced  by 
the  density  error  profile. 


a)  Noise  Error  in  Temperature 

The  error  in  temperature  due  to  the  noise  error  in 
density  is  given  by  Table  4,  line  1.  The  noise  error  in  temperature 
appears  in  the  computer  output  of  the  program.  The  unremoved 
bias  error  in  the  density  profile  due  to  smoothing  produces  a  bias 
in  temperature.  The  second  line  in  Table  4  shows  this  error. 

The  effect  of  a  ten  percent  error  in  the  initial 
temperature  on  succeeding  temperature  calculations  is  also  given 
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Table  4 .  Degree  Error  in  Temperature  for  Sphere  with 
Apogee  115  Km 


in  Table  4.  After  10  Km  of  flight,  the  effect  of  the  ten  percent 
error  nas  decreased  to  a  1°  error  in  temperature. 

The  effect  of  CD  errors  and  vertical  winds  on  temperature 
depends  upon  their  effect  on  pressure.  Since  the  bias  CD  error 
profile  and  vertical  wind  profile  are  unknowns,  the  resulting 
pressure  and  temperature  errors  cannot  be  estimated.  One  point 
of  interest  is  worth  mentioning.  If  all  drag  coefficients  are 
biased  by  a  constant  percent  error,  this  will  produce  the  same 
percent  error  in  density  and  pressure,  and  hence,  by  Equation  1, 
will  result  in  no  error  in  the  temperature.  It  is  a  changing 
bias  error  in  CQ  that  will  detrimentally  affect  temperature. 

Frequency  Response  of  Thermodynamic  Data 

The  amount  of  detail  that  can  be  observed  in  the  thermodynamic 
output  data  is  determined  by  the  frequency  response  of  the  density 
filter.  Table  5  gives  the  percent  of  the  amplitude  of  a  1,2,5, 

10  and  20  Km  sinusoidal  wave  in  density  that  is  retained  in  the 
output  data.  For  example,  at  70  Km  a  10  Km  density  perturbation 
with  amplitude  X  would  appear  in  the  computer  printout  as  having 
amplitude  .72X. 

Wind  Errors 

The  only  significant  sources  of  wind  error  are  the  noise  and 
bias  errors  in  velocity  and  acceleration  resulting  from  the  smoothing 
process.  For  a  sphere  apogee  of  115  Km,  the  one  sigma  noise 
error  in  a  wind  component  is  estimated  in  Table  6.  The  computer 
program  calculates  the  noise  error  in  each  wind  component  and 
includes  both  in  the  printout. 

Frequency  Response  of  Wind  Data 

The  bias  error  in  winds  introduced  by  the  smoothing  function 
depends  on  the  wind  field  itself  and  hence  varies  from  flight  to 
flight.  The  frequency  response  of  the  wind  filter  determines  the 
degree  of  bias  that  will  occur  for  various  wavelengths.  Table  7 
presents  this  frequency  response.  For  example,  from  Table  7  it 
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'•Alio'cr  APPUIUOt  CF  AKOOTHEO  0:K' I JV  W»  VE  '  1  tf  *  MPL  I  T  UOE  CF  CP  I  Gl  Nil 

WAVl  AS  »  FUNCTION  OF  AlTItuOE  ANC  VAVEICXITM. _ _ 

(/LTITLCE  ANC  WAVELENGTH  PEASLPEO  l«  X  I  LOPE  I  £  »S.  1 

^ALTlfuoJ 

"  "  X  IOC  X  90  «  ‘  PC  X  TO'  >  CO'  X  "  5C'  X  '  AC  'X  30  X  " 

XXXXXXXXXXXXXXX  xxxxxxxxxxxxx>>xxxxxxxxx»x»>xxxxxxxxxxxxxx 
X  X  X  X  X  X  I  X  "  X  . 

1  X  C.C1  X  C.01  X  O.Cl_X  J.0l_>  O.Cl  XC.22  X^.A2  X  0.75JC _ 

'x  X  X  '  X  X  'x  “  >  X  X 

2  X  C.CI  X  C.Ol  X  0.C1  X  O.OA  X  0.22  X  C.;2  X  O.BO  X  0.90  X  _ 

xxxxxxxxx 

5  X  C.CI  X  C.02  X  0.  C  3  X  0 .  A  7_  X  C.3A  X  C.90  X  C.96  X  0.99  X  _ 
...  X  X  '  X  "  X  x  x"'  "  x  '  X  X 

JO  X  C.A7  X  C.51  x  0.‘9  X  0.J2  X  0.97  X  0.9?  X  0.99  X  0.99  X _ 

X . X  X  x”  X  X  X  X  X' 

_  20  X  C.9A  X  C.9S  X  0.94  X_CK9  B_X_  C .  9  9  X_0.  S_9_  X_0.  99  J  0 . 9 ±  X  _ 


Table  5.  Frequency  Response  of 
Thermodynamic  Data 
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RATIO  Of  AXRLIUCE  Of_  S".C01fE0  Sr  l;3ID*L_WIN0  1C  tRPLITUOE  OF 
ORIGINAL  WINE  AS  A  I  USCT  1  ON  OF  Al  T  l  TUDc  ANC  I.  A  VELCNGTH. 

__ _ (XLTIRCE  ANC  WAVELENGTH  HE  »SUREO_  IN  KILOHcTERS.  I _ 

_ _ _ AIT  I  JUDE _ 

_ X..10C  X  90  X  fC  X  70  »_  tO  X  5C  X  _A0  X  30  _X 

xxxxxxxxxxxxxxx  xxxxxxxx'xxxxxxx  xxxxxxxxxx  x'xx  xxxxxxx  xxxxxxx 

A  _X  X  X  X  X  _X_  X  _  X  X 

A  i  X  O.CI  X  C.01  X  O.CI  X  0.01  »  0.01  X  O.C2'x'  C.  5  7  X  0  •  96  X 

y_x  x  x  x  x  x  i  x  x 

E  2  X  C.C1  X  C.01  X  O.Ci'x  0.03  X  0.10  X  0. A 0 *  X~  C." 76  X  0.99  X 

l _ X.  _  _  X_  .  X  X _  X  _X _  X X  X 

E  5  X  C.C2  X  C .  03  X  olct  X  0.11  li  OTlO  X  O'ZVj  X  0.99  X*  o.'99'x 

A  _  X  X  X  X  X  X  X  X  X 

t  10  X  C . 1 A  X  {.17  X  0.2A  X  0.85  X  0.97  X  0. 99  X  C.99  X  0.99'  X 

TX  X  X  X  X  X  X  X  X 

f  20  X  0.58  X  C.13  X  O.tt  X  0.90  X  0.9e  X  0.99  X  0.99  X  0.99  X 


Table  7.  Frequency  Response  of  Wind  Data 


is  observed  that  at  70  Km  a  10  Km  sinusoidal  wave  in  the 
field  will  appear  in  the  computed  output  to  have  only  85 
of  its  true  amplitude. 


wind 

percent 
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